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The two-point correlation function is calculated in the Lorentzian EPRL spinfoam model, and 
shown to match with the one in Regge calculus in a proper limit: large boundary spins j — > oo, and 
small Barbero-Immirzi parameter 7 — > 0, keeping the size of the quantum geometry A ~ finite 
and fixed. Compared to the Euclidean case, the definition of a Lorentzian boundary state involves 
a new feature: the notion of past- and future-pointing intertwiners. The semiclassical correlation 
function is obtained for a time-oriented semiclassical boundary state. 



I. INTRODUCTION 

Spinfoam amplitudes provide a covariant definition of 
the dynamics of Loop Quantum Gravity Q3-H]. A basic 
test the amplitude has to pass is that, its semiclassical 
behavior matches with classical gravity. The vertex am- 
plitude has been shown [5T-fT2] to correctly determine the 
Regge action for discrete gravity. Given this result, the 
next test for the theory regards the behavior of small 
quantum fluctuations around a classical solution. In the 
Euclidean theory, the two-point correlation function for 
the Penrose metric operator (also called the spinfoam 
propagator) restricted to a single vertex is calculated in 
[13] and the three-point function in [14J. In this paper 
we will consider the physically relevant Lorentzian case. 

In this paper we compute the Lorentzian two-point cor- 
relation functions for the Penrose metric operator. The 
setting is the one introduced in [15] and developed in 
13, 16 24 . In particular, we restrict to a single spin- 
foam vertex and compute correlations on a semiclassical 
state peaked on the spacelike boundary geometry of a 
Lorentzian 4-simplex. For a discussion of the role of this 
truncation see |2"5] . 

Here we briefly review the procedure for completeness. 
We consider a 4-dimensional manifold 1Z homeomorphic 
to a 4-ball, with its boundary given by a 3-dimensional 
manifold £ = dlZ homeomorphic to the 3-sphere S 3 . We 
associate to S a boundary Hubert space of states: the 
LQG Hilbert space %•£ spanned by (abstract) spin net- 
works. We call a generic state in H^. A spin foam 
model for the region 1Z provides a map from the bound- 
ary Hilbert space to C. We call this map (W\. It provides 
a sum over the bulk geometries with a weight that defines 
our model for quantum gravity. The dynamical expecta- 
tion value of an operator O on the state |^) is defined 
via the following expression 



This expression corresponds to the standard definition in 
(perturbative) quantum field theory where the vacuum 
expectation value of a product of local observables is de- 
fined as 



(O( Xl )---O(x n )) 
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The vacuum state Vf'o^] codes the boundary conditions 
at infinity. The operator O can be a geometric operator 
as the area, the volume or the length [2"6"rf3"3"] . The geo- 
metric operator we are here interested in is the (density- 
two inverse-) metric operator q ab (x) = 5 13 ' Ef{x)E b (x). 
We focus on the connected two-point correlation func- 
tion G abcd (x,y) on a semiclassical boundary state l^o)- 
It is defined as 

G abcd (x, y) = (q ab (x) q cd (y)) - (q ab (x)) (q cd (y)) . (3) 

We generalize the semiclassical states used in |13j . 
which is obtained via a superposition over spins of spin- 
network states having nodes labeled by Livine-Speziale 
coherent intertwiners [3~4Tj56"|, with coefficients proposed 
by Rovelli in [TS] . The origins of this semiclassical state 
is studied in |37j . Here we generalize this semiclassical 
state to be Lorentzian, in the sense that it is related to 
Lorentzian geometry. 

We assume that the Barbero-Immirzi parameter 7 is 
positiveQand restrict ourselves to the first order of vertex 
expansion, i.e., we consider the spinfoam with a single 
vertex. Our main result is the following. We consider the 
limit, introduced in |13) and discussed in [51 IIP), where 
the Barbero-Immirzi parameter is taken to zero 7 — > 0, 
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and the spin of the boundary state is taken to infinity 
j — > oo, keeping the size of the quantum geometry A ~ 
jj finite and fixed. This limit corresponds to neglecting 
Planck scale discreteness and twisting effects, at large 
finite distances. In this limit, the two-point function we 
obtain exactly matches the one obtained from Lorentzian 
Regge calculus 38 . We therefore extend to Lorentzian 
signature the results of [15] . 



II. LORENTZIAN GEOMETRY AND 
LORENTZIAN SEMICLASSICAL BOUNDARY 
STATES 

Semiclassical boundary states are a key ingredient in 
the definition of boundary amplitudes. Here we describe 
in detail the construction of a boundary state peaked on 
the intrinsic and the extrinsic geometry of the bound- 
ary of a Lorentzian 4-simplex. The construction follows 
the Euclidean case |13) , by considering Lorentzian geom- 
etry of a 4-simplex with space-like tetrahedra. As in the 
Euclidean case |13j . it uses the Livine-Speziale coherent 
intertwiners [33] (see also [35]) together with a superpo- 
sition over spins as done in jTS], [TC] , and thus it can be 
considered as a Lorentzian version of an improvement of 
the boundary state used in [2"0"H2"2"] where Rovelli-Speziale 
gaussian states [39] for intertwiners were used. 

We consider a simplicial decomposition A5 of S 3 . The 
decomposition A5 is homeomorphic to the boundary of 
a Lorentzian 4-simplex: it consists of five tetrahedra a 
which meet at ten triangles (ab) (a,b = 1...5 and a < 
b). In particular ,we focus on the Lorentzian geometry, 
where all the five tetrahedra a are space-like, i.e., all the 
normals to the tetrahedra are time-like. Suppose the 
time- like normals to the tetrahedra are outward-pointing, 
for convenience. Then for each triangle (ab), there is a 
corresponding "wedge" composed by the two tetrahedra 
a and b which meet at the triangle (ab). Since all the 
tetrahedra are space-like, the tetrahedron comes in two 
types: the outward normals are either future-pointing or 
past-pointing. The wedges are then classified into two 
types: it is called in [7J 130] thick wedge if the incident 
tetrahedra are of same pointing type, which means both 
future-pointing or both past-pointing, otherwise called 
thin wedge. We introduce a quantity H a b to denote the 
Lorentzian geometry of the triangle (ab) in the following 
way: 



J 0, thick wedge 
[ 7r, thin wedge 



Now we consider the complete graph T§ dual to the 



decomposition A5: 



V5 




t'2 



with five nodes v a dual to tetrahedra of the 4-simplex and 
ten links l a b, (a < b) dual to the corresponding meeting 
triangles. Consider the group SU(2) and associate an 
irreducible representation ji (a spin) to each link I of the 
graph r, and an SU(2) intertwiner i n to each node n of 
the graph. The spin network states \Ts, j a b, ia) supported 
on this graph span the truncated SU(2) Hilbert space 
Ur 5 of LQG. 

Now let us introduce a (overcomplete) basis of 
Hr 5 , which is the coherent spin network states with nodes 
labeled by Livine-Speziale coherent intertwiners. We 
start from associating an SU(2) coherent state \j, fi(£)) 
[3j] [31] in the spin-j representation to each link of the 
graph T5, labeled by the spin-j and a unit vector n(£) 
defining a direction on the sphere S 2 , associated to a 
normalized spinor ( e C 2 . In fact, one can obtain these 
coherent states from the maximal weight vectors \j,j), 
when m = j on the basis states \ j, to), which minimize the 
(SU(2) invariant) uncertainty A = \(,P) - (J) 2 \ in the 
direction of J z . Starting from the highest weight, an in- 
finite set of coherent states on the sphere SU(2)/U(1) ~ 
S are constructed through the group action, \j, n) = 
n li>j)) where n is a unit vector defining a direction on 
the sphere S 2 and n an SU(2) group element rotating 
the direction z = (0, 0, 1) into the direction ft. Just as 
I j, j) has direction z with minimal uncertainty, \j, ft) has 
direction ft with minimal uncertainty. One can go fur- 
ther to use a normalized spinor £ g C 2 to label an SU(2) 
group element 

-<«=(! t) 

and the corresponding vector ft(£). 

These coherent states are peaked on the geometry of 
a classical triangle: n are associated to unit-normals to 
triangles of a tetrahedron, and j areas of the triangles. 
We assume all the normals outward to the tetrahedron, 
which satisfy 

jl«4 + J2«2 + J3«3 + jA = 0, (6) 

thus we associate to each triangle (ab) normals — n a b and 
Hba when a is target of the triangle and respectively b is 
the source. The antipodal vector — n(£) can be associated 
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to J£, i.e. — H(£) — n(J£), with 



J 



Calling t;™ 1 " ™ 14 the standard recoupling basis for inter- 
twiners, we can define the coefficients 



(7) 



Ti(ni, 



-* \ mi ■■■m,A '\r /-* 

.,n i ) = v i 4 T mi ... m4 (n 1; 



.,n 4 ). (14) 



A Euclidean coherent intertwiner between the repre- 
sentations ji, . . , j'4 is defined as 



We define a Lorentzian coherent spin network 
1^5, jab, T a ) as the state labeled by ten spins j i, 
and 4x5 normals n j and given by the superposition 



r 4 5 

$(n u ..,n 4 ) = j ^ 2 d/i n0'a,m a |%a,n o ) (8) |r 5 , j a6 , T a (n)) = £ ( J] T ia (n a6 )) |r 5 ,j a6 , i a ) 



up to a normalization, with the SU(2) coherent states 
\j, ft) satisfying the closure constraint ([6]). Here we use h 
to denote the SU(2) group element as well as the corre- 
sponding representation. Let us now introduce the notion 
of Lorentzian Livine-Speziale intertwiner, which general- 
izes the Euclidean Livine-Speziale intertwiner j34j to the 
one related to the Lorentzian geometry. Since one has the 
closure condition for the 4-simplex and all the tetrahedra 
are space-like, with outward normals, there is at least 
one (and at most four) tetrahedron with past-pointing 
outward normal. The Lorentzian coherent intertwiner 
T(ni, . . ,714,) can be obtained from the Euclidean one Q 
by time reversing the normals of past-pointing tetrahe- 
dra: 



• For future-pointing tetrahedra, 

T(ni, . . ,n 4 ) = $(ni,. . ,n 4 ); 



(9) 



• For past-pointing tetrahedra, 

r(n 1 ,..,n 4 )=T$(-n 1 ,..,-n 4 ). (10) 

The effect of time reversal T on the coherent states are 
given by 

T\j,n) = (-iy\j,-n), (11) 
and thus on the (past-pointing) coherent intertwiner by 



T$ (7H,..,n 4 ) = (-lJ^W-^Hfi, 



-n 4 ) 



(12) 



Here we assume the normal to the 5th tetrahedron are 
future-pointing and thus a = 1, ...4. The minus signs 
in the r.h.s. of equation ( 10 ) keeps normals outward- 
pointing. Combining the two components ^ and (10 1, 
one obtains 



Y Q (ni, . . ,n 4 ) =exp(-i^n a6 j ab )$ a (ni,..,n 4 ), (13) 

b>a 

with H a b related to Lorentzian geometry in equation 
Q. From the intertwiner relation (13), one can 
call exp (— i ^ a<f) n a fcj a f,) Lorentzian- geometry phase, 
which maps the Euclidean coherent intertwiners to the 
Lorentzian ones. 

This Lorentzian coherent intertwiner can be used to 
label nodes of Lorentzian coherent spin network state. 



(15) 

If one uses \T 5 ,j ao , $ Q (n)) to denote the Euclidean coher- 
ent spin network state with nodes labeled by Euclidean 
coherent intertwiner $ a (n), as in |13) . the Lorentzian 
coherent spin network state jTs, j a b, Y (n)) can be ex- 
pressed as |r 5 , j a b, $ a (n)) times a Lorentzian-geometry 
phase: 



|r 5 , jab, ^ a {n)) = cxp M^n Qb j Qb )|r 5 , j ab , $„(*?)). 



a<b 



(16) 



The two phases from the intertwiner relation (13) and 



from the spin network relation (16 
sense that a is fixed in the former 



are different, in the 
:>ut free in the latter. 



In fact, the latter phase is product of the five phase in 
the former form. However, here, we call both of them 
Lorentzian-geometry phase, if there is no confusion rises. 

A Lorentzian semiclassical state peaked both on intrin- 
sic and extrinsic geometry can be given by a superposi- 
tion of Lorentzian coherent spin network states: 



l*o>=$>i..*»C?')|j.T(rr)) 



(17) 



with coefficients i/jj o 
phase, 



,(j) gi ven by a gaussian times a 



VWoC?) = eX P (- i ^2^f ( t ) o b (jab ~ (jo)ab) I x - 
ab 

(ab)(cd) 3 oh - (Jo) - (jo. 



x exp 



ab.cd 



cd 



y/(jo)ab \/ (jo)cd 



(18) 

where 4> labels the simplicial extrinsic curvature, which 
is an angle associated to the triangle shared by the tetra- 
hedra; the 10 x 10 matrix a^ ab ^ cd ^ is assumed to be com- 
plex with positive definite real part. These coefficients 
are originally proposed in [15] and used to calculate the 
Euclidean two-point [13] and three-point [14] correlation 
functions. Here we will use the Lorentzian semiclassical 



state (17) to calculate the Lorentzian two-point correla- 



tion function. 



III. THE EPRL AMPLITUDE IN LORENTZIAN 
THEORY 

In this section, we give a brief introduction to the 
SL(2, C) EPRL amplitude of a coherent spin net- 
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work. Throughout this paper, SL(2, C) refers to the 6- 
dimensional real Lie group of 2 x 2 complex matrices 
with unit determinant, and is called simply the Lorentz 
group. It covers the group of proper orthochronous 
Lorentz transformations, SO + (3,1), which is the com- 
ponent of the group 0(3, 1) connected to the identity. 

The principal series of irreducible unitary representa- 
tions of the Lorentz group SL(2, C) are labeled by two 
parameters (k,p), with k an integer and p a real num- 
ber |42) . Given a carrier space the canonical ba- 
sis is given by the basis diagonalizing simultaneously the 
Casimir operators, which is denoted as \(k,p); j,m). 

The SL(2, C) EPRL amplitude of a single 4-simplex for 
a Lorentzian coherent spin network state \T^,j a b, T a (n)) 
reads 

(W\T 5 ,j ab ,T a (H)) 

= e -ii: ab n abjab J Y[dg a l[P ab (g), (19) 



SL(2X) 5 



a (ab) 



with 



Pab(g) = (jab,-Kam\y^9a l 9bY\jah,riha{0)- (20) 

Notation is as follows. The indices a, b = 1, 5 label the 
tetrahedra on the boundary of the 4-simplex and (ab) la- 
bels the triangles between the corresponding tetrahedra; 
the integral is over one group element of SL(2,C) per 
each tetrahedron. We restrict ourselves to the spacelike 
tetrahedra. We use g to denote the group elements, as 
well as the corresponding representations. The EPRL 
embedding map Y embeds the spin-j irreducible repre- 
sentation Hj of SU(2) to the irreducible unitary repre- 
sentation H(k,p) °f SL(2,C), given by 



Y\j,m) = \(j,jj);j,m). 



(21) 



Note that we associate to each triangle (ab) normals — n a b 
and Hba when a is target of the triangle and respectively b 
is the source. That is why we have a minus sign in the bra 



coherent state in the definition of P ab in equation (20 1 
And the Hermitian inner product is defined as (z, w) = 

Z W + Z-lW-l. 



To see this amplitude ( 19 ) explicitly, one can turn to 
a representation of the Lorentz group SL(2, C) on the 
space H(k,p) °f homogeneous functions of the complex 
affine plane C 2 — {0,0}, 



f(az)^ = a- 1+vp+k d- 1+lp - k f(z) 
Va G C - {0}, 



(fc,p) 



with the group transformation 

g: f(z)^f(g T z). 



(22) 



(23) 



The canonical basis is denoted as /^(z)( fe ' p ). And the 
inner product is given by 



(f,g) = J dz fg (24) 

with dz = |(zodzi — zidzo) A (zodzi — Zidzo). This inte- 
gral is invariant under the scaling z — > az, according to 
the homogeneity ([22]). To modulo this equivalence rela- 
tion, one can choose ip(z) = f(z, 1) associated with each 
f(z ,zi) G %(k,p)] the functions <f(z) forms a realization 
of T~L(k,p)> which we can still call H(k jP )- Functions ip(z) 
can be considered as the homogeneous functions on the 
complex projective line V, the subspace of the complex 
affine plane, modulo the equivalence relation az = z. For 
calculating simplicity, we will still keep the formulae of 
f(z) on the complex affine plane in the following and 
reduce to f(z) G V when necessary. 

In this representation of homogeneous functions, 
SL(2,C) coherent state \(k,p); k, n(£)} with lowest spin 
k can be written as [7] 



\ ip—X—h 



2A 



(25) 



And hence equation ( 20 1 can be rewritten [43] as 



Pab =(jab, ~n a b(0\ yt 9a 1 9b Y\jab, ™ba(£)) 

= ((jabn3ab); Jab, -n ab (0 \ 9a 1 9b \(jab,13ab)\ jabrfbaiO) 

= [ dz g a ff \(z)(^ra ab ) gb /£* ( z )(f.»,7i.6) 
Jcr 1 

= J dz ffl b (9 T a z ) Uab '' 10ab) ft (9^) iub ^ b) 
77 J CP 1 

= ^ f dz fe^)"^^-^^^ 
n Jcr 1 

_ dj ab f f (Zba, Zba)\ 1 ~ a " b f (J£,ab, Z a b) 2 (Zba,{;ba) 2 \'' ab 

K Jcr 1 \(Z a b, Z a b) J V (Zab, Z a b)(Zba, Zba) ) 
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where dz ab = -((Z ab , Z ab )(Z ba , Z ba )y 1 dz, Z ab = g\z 
and Z ba = g' b z; (, ba and J^ ab are spinors associated re- 
spectively with ri ba (£) and — n ab (£), as introduced in 
equation ([7]); note that g is used to denote the group 
elements, as well as the corresponding unitary represen- 
tations; the property of unitary representation is consid- 
ered in the 3rd step , equation ( |25[ ) used in the 5th step 
and in the 6th step, the integral variable z is changed 
into its complex conjugate z, where the minus sign comes 
from the integral measure. If we denote the integrand in 
equation (26) as K ab : 



I - i I (Zba, Z ba ) \ 7:,ab 

K ab {g,z) = ( ^-r I 

\Aj6i ^ab) , 



( (J£ab> Z ab } 2 (Z ba ,^ ba ) 
V (Zab, Z ab )(Z ba , Z ba ) 



2\ 3ab 



the integral (26 1 is simply expressed as 



Pab = 



'-dz ab K ab 



(27) 



(28) 



Thus the EPRL amplitude ( 19 ) can be written as 
(W\T 5 J ab ,T a (n)) 

a—l \a<b / 



(29) 



where the "action" S is given by 

(J£ab, Z a b) 2 (Z b a,£,ba} 2 

y Jab iUJd, 

a<b 



-ijj ab lOf 



\Z a b, Z ab )(Z ba , Z ba 
(Z ba , Z ba X 



(Z a b, Z a b) 



iRabjab 



(30) 



Note that the "action" S here has one more term 
— i^2Habjab than the one in [7], since we consider the 



Lorentzian coherent spin network states. This expres- 
sion, however, is ill defined, due to the fact that the inte- 
gral may diverge. This issue is addressed and answered 
in [33J 05] , where it is shown that the source of the diver- 
gence is a redundant integral over SL(2, C) in the vertex 
amplitude ( [29] ) . It is then immediate to regularize the 
vertex amplitude by removing one SL(2, C) integration. 
The resulting amplitude with an integral over Ila=i d<?a 
is proven in [H, 45J to be finite. 



IV. LORENTZIAN TWO-POINT FUNCTION 
AND ITS INTEGRAL FORMULA 

Following [T3], the connected two-point correlation 
function G^f on a semiclassical boundary state |^ ) 
is defined as 

rtabed / ™ jpb rpc rpd \ / rpa Tpb\ / jpc rpd \ /qi \ 

G nm - {^n^n E m -E m ) ~ K^n ^n) \ fc m' ' &m) i (31) 

where (E%)i is a flux operator through a surface f an dual 
to the triangle between the tetrahedra a and n, parallel 
transported in the tetrahedron n. Here the dynamical 
expectation value of an operator O on the semiclassical 
state I*,,) is defined via 



(O) 



(gw 

(W\9 ) 



(32) 



The semiclassical state \& )> which is introduced in the 
end of section [TTJ can be simply expressed as a superpo- 
sition: 



l*o>=5>^'> T (")>> 



(33) 



with ipj given by equation ( 18 ). Thus the two-point func- 



tion (31 ) can be also written as a superposition: 



G 



abed 



E^ 3 (W\K-K E^Ef n \j,T(n)) E^jiWlK-E^Tin)) Y. ] ^(W\E< m -E d m \ ]l T(n)> 



T,i^(w\j,r(n)) 



E^j(W\j,T(n)) 



(34) 



To see this explicitly, let us go first to derive integral 
expressions for (W\E« ■ E b n \j, T(n)) and (W\E% -E h n E c m - 
E^lj, T(n)). As in [T3], one can introduce an "insertion" 



and obtain the integral expression of (W|.E^ -E^\j, T(n)) 
and (W\E%-E* E^-E^j, T(n)) in terms of this insertion: 



lab 



(jab,~n ab (0\Y^g- 1 g b Y(E-y\j ab ,n ba (0), (35) 
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(W\E a n -E h n \j,T{n)) = 



j ex P (-< ]r n abJab ) n d 5a ^QL^ b n' pcd 

ab a—1 cd 

(W\E«-E h n E c m -Ei\j, r{n)) = / cxp n o6 j a6 ) J] d. 9a fyQ^Q^wQLQJmJI P ' 



(36) 



cd 



where the product Y[ is over couples (cd) different from 
(na), (n6), (mc), (md). 



Now let us come to express the insertion Q\ in ( 35 1 as 



map Y 



y J l ab\jab,m ab ) = Ja b Y\jab,m ab ) 



(37) 



and the fact that the generator J l ab of SU(2) can be ob- 
tained as the derivative 



, d 



a l =0 



ab± , 



a group integral. Using the invariance properties of the we have that 



(38) 



where 



Using 



T ab ={j ab , -n ab (0\ Y^g^gbYiE^lU, n ba (0) 
. d 



da 1 

d_ 
da 1 



a i=° Jcp 1 



dz g a fjfjz)^^ g b e- M ' r - fi*( z )V*™*) 



-dz. 



(Zba, Z, 



ba 



a I=0 y cp l 7T (Z ab ,Z ab )(Z ba ,Z ba ) \{Zab,Z ab ) 

I 



({Jiab, Z ab ) 2 (Z ba , £ ba ) 2 \ 
(Z a b, Z ab )(Z ba , Z ba ) J 



equation ( 39 ) turns out to be 



Z ba = (9b e~ iam )h ab = e 



—lOiiTi '/ 

Aba 



{Zba, Zba) ~ (Zba, Zba) 



(40) 



(41) 



(39) 



Let 



=«7 



n Jcp 1 ah \{Z a b,Z a b) 

djab f d ~ / (Z b g, Z bo } X 
f JCP 1 \(Z a b,Z a b 
2jqfe 

(Zba,£,ba) da 1 



a 

da 1 



((J£ab,Z a b) 2 (e laiT Zba,£,ba) 2 \ 
(Z a b, Z ab )(Z ba , Z ba ) J 

(Jtjab, Z ab ) 2 (Z b , t; ba ) 2 \ 3ab 
(Z a b, Z a b)(Zba, Z ba ) 

Z ba , £,ba) 



_ d Ub 



dz 



ab 



( e la ' T Z ba ,iba) 

f (z bai z ba ) \ i ^ b f (Ju,z a b) 2 (Zb,iba) 2 \ 3ah b <y 

\(Z a b, Z a b) ) \ (Zab, Z a b) (Zba, Zba) ) (-2 



Zba,£,ba) 
(Zba, £ba) 



(42) 



Kb = Uab 



(Zba,C 



,ba I 



(43) 



and use -K" (, given by equation ( 27 1 , one can rewrite equa- 
tion ( 42 ) simply as 



Tab = / ^d^KaUAahT. 



(44) 
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This equation can be used to rewrite equations ( 36 1 as 



(W\E^.E b n \j,T(n))= [Ud 9a [[Yl %^dz Q ^) qf e s 

J a=l J \a'<b' 77 / 

(W\E« ■ E b n E? n ■ Ei\T(n)) = j ' f[dg a j (\[ ^dz Q , b , j qfqf. 



(45) 
(46) 



with 



ab _ a a Ab 

Hn — n r; 



I 

Here we remove a redundant 5g§ integral, as discussed 
in Sec. Ill Then the two-point function (34) can be 



(47) reexpresscd in terms of group integrals: 



abcd Ej i>J I ^9 S w z q a n b q c Je s Ej i>J I <5 4 ff S w z q a n b e s £ ^ J 5±g S w z q% 



s~iabcd 



(48) 



where the group integral is over 5 g = Ila=i an d S 10 z thus the two-point function (48) can be reexpressed as 
is short for the integral measure n„<& ^dzah over V. G ( A )-. To study large-spin limit turns then to study large- 

A limit, via stationary phase approximation. In Sec. |V A| 
we give a brief framework of this technique. Then in Sec. 
V B||V C we give the detailed calculation. 



V. ASYMPTOTIC EXPANSION OF THE 
TWO-POINT FUNCTION FOR LARGE SPIN 

In this section we study the large- j Q asympototics of 
the correlation function ( 48 1 . We use the technique devel- 
oped in [T3]. The idea is rescaling the spins j ab and {j )ab 
by an integer A so that j ab -> Xj ab and (j ) a b -> X(jo)ab, 



A. The rescaled correlation function and 
stationary phase approximation 

As in [T3], let j ab -> Xj ab and (j ) ab -> X(j ) ab , we 
rescale the correlation function ( 48 ) as 



rrtcdr^ EjJ d^gS^zqfq^e^ £. J d^gS^zqfe^ ^ I <fig8™z<tfe**>* 



Ej/ d i gS 10 ze xs ^ Ej / d 4 g8 w ze xs ™ EjJ d i g8 w ze 



ASto 



(49) 



where the "total action" is defined as Stot = log ip + S or 
more explicitly as 



Stot(j)5) z ) 



2 ^ 

ab,cd 



(ab)(cd) Jab {jo)ab jed (jo) erf _ . , a b 



cd 



i }^ 70o {jab ~ UoU) + S(j, g, z) (50) 

ab 



Using Euler-Maclaurin formula, one can evaluate the 
sums over spins j using integrals in the large A limit: 



d 10 J qfe xs ^+0(\- N ) 

> , (51) 
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so that the rescaled correlation function (49) can be 
approximately expressed in the large A limit 



.,„,„,, JS 10 j d A g5 10 zq% b q%e xs ™ J S w j d 4 g S w z qfe xs ™ J S w j d A gS w zq c Je XSi - 



Gabcd ( \ \ 
nm l A J 



JS 10 j d i g5 10 ze xs f 



J6 w j d 4 gS w ze xs ^ J S w j d i g5 w ze xs ^ 



(52) 



r 



We will study the large-A asymptotics of expression (52 1. Let us rewrite the two-point function (52 1 formally as 



G(\) = 



J dx p(x)q(x) e xs( - x ^ J dx p(x) e xs ^ J dx q(x) 



Jdx 



fdxe xs W ' fdxe xs W 



(53) 



then the asymptotic expansion of G(A) for large A is given 
by 



1 



G{\) = -{H- 1 r P ' i {x )q'Ax ) 



A 



(54) 



Here x is the critical point, i.e. the stationary point 
where the real part of the action vanishes, KeS(x ) = 0; 
p\ = dp jdx 1 , H is the Hessian matrix at the critical point 
H = S"(x ). Our task is to obtain the critical point, the 
derivative of the insertions in Sec. IV Bl and the Hessian 



in Sec. WU 



B. The critical point and the derivative of 
insertions 



The critical point is the one where the real part and 
the derivatives of the total action 5 tot vanish. The real 



part of the total action ( 50 ) is given by 



ReS to t = - 1] 7(Rea)( a 



b)(cd) Jab — (jo)ab jcd — (jo)cd . , | (J£,ab , Z a b) | 2 | (Zba , £ba) \ 2 
r ^ab log 



ab,cd 



\J (j'o)afe V {jo)cd 



(ab) 



(Z a b, Z ab ) (Zba, Zba) 



(55) 



r 



Having assumed that the matrix a in the boundary state 
has positive definite real part, we have that the real part 
of the total action is negative or vanishing, ReStot < 0. In 
particular the total action vanishes for the configuration 
of spins jab and group elements g satisfying 



jab (jo)ab , 



J£ab 



Zab 



Jab, 



-,ba 



Zba 



J ba, 



(56a) 



(56b) 



where || Z a b || is the norm of Z a b induced by the Hermi- 
tian inner product, and (f> a b and (f>b a are phases. 

The requirement that the variations of the total action 
with respect to the spinors z a b and z a b vanishes, ^z^tot = 
SzS to t = 0, lead both to 



,g<iJ£,ab 



gb£,ba 
\\Zba\\ 



(57) 



evaluated at the maximum point (56b). For the group 
variables, S g Stot = leads to 



^ jabn ab = 

b:b^a 



(58) 



evaluated at the maximum point ( 56b ) . In fact the nor- 



mals n a b in the boundary state are chosen to satisfy the 
closure condition at each node. Therefore the critical 
points in the group variables are given by all the solu- 
tions of equation (56b). 



The variations of the total action with respect to the 
spins j turns out to be 



dS t , 



(ab)(cd) 



djab 



c d y/ {jo )ab \J {jo ) cd 



{jcd {jojed) _ ^0tt6 



OS 

djab ' 

(59) 
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Imposing the maximal-point equation ( 56 1 and the 



critical-point equation (571, equation (59) is reduced into 



djab 



(60) 



where the parameter /i — ±1 measures the discrepancy 
of the orientations of the 4-simplex a: there are two 



orientations of the 4-simplex er, one inherited from the 
Minkowski space where a is embedded, the other induced 
from the boundary data; fi = 1 if these two agree and 
fj, = — 1 otherwise. The requirement that equation (60) 



vanishes selects fi = 1, which means we only consider the 
case when the orientation of the boundary data agrees 
with the one induced from the Minkowski space. 

We end this subsection by the first derivative of the 
insertion qn b (g,z) evaluated at the critical point: 



<rr ab 



crit 



"T (jo)na (jo) n h 



\\Zar 



Han ' ^bn 



=l 2 {jo)na{jo)nb (j4n(£<Wi! LcT^n) - -(n an ' n bn )n, 
=J 2 {jo)na(jo)nb ( «L (£anK<7i£ a n) - ^ ("an ' n bn )ri a 



(Jo)ab 



(61) 
(62) 

(63) 
(64) 
(65) 



Hessian matrix of the total action 



(10 + 24 + 20) x (10 + 24 + 20) matrix 



Following the stationary phase approximation intro- 
duced in section IV Al once the Hessian matrix is ob- 



tained, one can get asymptotic expansion of the two- 
point function (|48|) by using equations (l54|) and ([61) 



Now let us come to calculate the Hessian matrix. 

The Hessian is defined as the matrix of the second 
derivatives of the total action where the variable g§ has 
been gauge fixed to the identity. We split the Hessian 
matrix into derivatives w.r.t. the spins j, w.r.t. the group 
elements g and w.r.t. z. The Hessian will then be a 



as 



"tot — 



SiS„S. 



024x10 
020x10 



j u g°s to 



hox24 0iqx2O 



H 9 



Sj5 z S St , 



0. 



(66) 



(67) 



We will now describe the non-vanishing blocks of this 
matrix. Qjj is a 10 x 10 matrix containing only deriva- 
tives with respect to the spins j ab , with elements 



Q(ab){cd) — Sj ab Sj ci S t ot 



(ab) (cd) 



crit y/(jo)aby/(jo)cd 



S- 



crit 



Rcggc ■ 



(68) 



r 



H gg is a (4 x 6) x (4 x 6) matrix containing only deriva- 
tives with respect to the group elements g a . Note that 
due to the form of the action, derivatives with respect to 
two different group variables will be zero and it will be 
block diagonal 



performed by splitting the SL(2, C) element into a boost 
and a rotation generator. This gives 



H 99 = 

















#22 














#33 














H44 



(69) 



Each H aa is a 6 x 6 matrix. The variation has been 



rrrr Tjbr 
U ( al )(aj) U (ai)(aj) 

n {ai)(aj) U (ai)(aj) 



(70) 
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with 3x3 matrices 



spect to the spinors z ab and z a b 



H lii)w) =\ E + nV + ieij ^ k ) ( 71 ) 



Tjrb 

U {ai){aj) 



H, 



(ai)(aj) 



rrbb 

n ( al ){a 3 ) 



i 



J2 3ab{S ij + rtv? + ie^ k n k ) (72) 



b:b^a 



- g E J'«*(-<^' + n ' ni + ie *W) (73) 

b:b^a 

i 
2 



=2 £ 3ab (l + 5 l] + n'n? + ie« \n k 



b:b^a 



(74) 



H, 



q q 

zz u zz 

n"_ a"__ 
zz zz 



(75) 



Since each spinor z ab has two compononts (z zi), each 
block S of H zz seems to be 20 x 20 matrix and thus 
H zz seems to be 40 x 40; however, this 40 x 40 matrix is 
degenerate due to the homogeneity of the representation 
functions ( 22 ) which we discuss in section III Remind 
that although we still keep the formulae of /(z), we have 
chosen a section for (z a z\) as (z /z\ 1). Hence the 
matrix H zz reduced to be 20 x 20 and non-degenerate, 
by removing the derivatives with respect to the second 
component z\. 



H zz = 



0i 



Oio 



xlO 



(76) 



H zz is a matrix containing only derivatives with re- 



with S zz and S zz vanishing, S zz and S zz diagonal matri- 
ces. The diagonal elements are first components of the 
matrix 



7ab( 



^{9aJ£,ab){g a J£,abY + («7 ~ l)ffaffl (*7 + l)5&Sb 



'Jab 



\z, 



bit 



(77) 



H gz — [Sgz Sgz) is a 24 x 20 matrix, containing 
derivatives with respect to spinor z ab and group element 
g a - Again, we consider the reduced matrix, by removing 



the derivatives with respect to the sencond component of 
spinor z a b- The non- vanish elements are the first compo- 
nents of 



H&Kab) = e#ab p^ ((! - H)^9a JUV - i(l + l)n a b(9a J U V) 



D. The 7^0 limit 

Now we obtain the Hessian matrix. This matrix is 
very complicated and its inverse is not easy to express 
explicitly. Despite this difficulty, it is possible to look 
at the approximate inverse in a proper limit, where the 
Barbero-Immirzi parameter is taken to zero 7 — > 0, and 
the spin of the boundary state is taken to infinity j 00, 



I 

keeping the size of the quantum geometry A ~ jj finite 
and fixed. In this limit, the Euclidean two-point [12] as 
well as three-point [13] correlation functions matches the 
ones in Regge calculus [Hj. More recently, Magliaro and 
Perini propose in [10] to understand this limit as the 
continuous spectrum limit of the fundamental geometric 
operators, since the Barbero-Immirzi parameter 7 enters 
the discrete spectra of area and volume operators |26l - 



11 



Now let us come to study the Hessian matrix 
in this limit and the correlation function as well. 

Before we take this limit, let us comment more on 
j — ^ oo. In the Lorentzian signature, we do not have 
the regular 4-simplex, so that we do not have j a b — j or 
(jo)ab — jo for all triangles. Since all j a b — > oo in the 
large- j limit, however, we can factorize each spin into 
jab = j^ab and let j — > oo and e a b finite. Thus we have 
a fundamental area A = jj and e a b is obtained when we 
use the fundamental area A — jj to measure triangle 

Now let us take the limit for the Hessian matrix. 

• The derivatives with respect to group elements 
H gg = jHgg, with 24 x 24 matrix Hg g independent 
from j and 7. 

• The derivatives with respect to spinors H zz = 
j(H e zz + 0(7)), with 20 x 20 matrix H% x indepen- 
dent from j and 7. 

• The derivatives with respect to group elements and 
spinors H gz = j(H gz + 0(7)), with 24 x 20 matrix 
Hg Z independent from j and 7. 

Thus we can express the Hessian matrix in the limit as 
Qu 



j] "10x44 
044X10 j(ff £ + 0(7)) 



(78) 



with 44 x 44 matrix H e independent from j and 7. The 
inverse of this Hessian matrix ( 78 ) is simply 



0l0x44 
;-l/YZTe\-l 



044x10 r l ((ff e )- l + o( 7 )) 



(79) 



Substituting (61)-(65) and the approximate hessian 



matrix (|78| into (|54|) we obtain the asymptotic expan- 

dqf dq% , 



sion of the correlation function 

fiabcd \ ^ r^ — 1 



p,q,r,s 

„4„-3 



(80) 



with X e independent from the boundary spin j and from 
7, as well as the combination 



R 



abed 
nm 



l 



7 3 Jo 



E 



dq: 



ab 



,<■<{ 



)<g,r<s 



' {pq){rs) dj pq dj 



(81) 



In terms of these quantities we have the following two- 
point correlation function 



= (7Jo) + 0( 7 )) 



(82) 



If we take the classical limit, introduced in [13] . where 
the Barbero-Immirzi parameter is taken to zero 7 — > 0, 
and the spin of the boundary state is taken to infinity 
j — > 00, keeping the size of the quantum geometry A ~ jj 



finite and fixed |^] the two-point function ( 82 ) we obtain 
exactly matches the one obtained from Lorentzian Regge 
calculus [55] , 

Similarly, higher order rt-point correlation functions 
can be computed using the generating functional tech- 
nique discussed in [49J [50] . 



VI. CONCLUSION 

In this paper we have studied in the Lorentzian signa- 
ture the two-point correlation function of metric opera- 
tors of EPRL spinfoam model. 

The analysis presented involves two distinct ingredi- 
ents. The first is a setting for defining correlation func- 
tions. The setting is the boundary amplitude formal- 
ism. It involves a boundary semiclassical state which 
identifies the regime of interest, loop quantum gravity op- 
erators E® ■ E h n which probe the quantum geometry on 
the boundary, a spin foam model (W\ which implements 
the dynamics. The formalism allows to define semiclas- 
sical correlation functions in a background-independent 
context. The second ingredient consists in an approxi- 
mation scheme applied to the quantity defined above. It 
involves a vertex expansion and a large spin expansion. 
It allows to estimate the correlation functions explicitly. 
The explicit result can then be compared to the correla- 
tion function of Regge calculus. In this paper we focused 
on the lowest order in the vertex expansion and the lead- 
ing order in the large spin expansion. In the following we 
collect some remarks on these two ingredients. 

The boundary semiclassical state |^ ) used here is 
Lorentzian, in the sense that it is related to Lorentzian 
geometry of 4-simplex. It is obtained via a superpo- 
sition over spins of Lorentzian coherent spin network 
1^5, jab, ^ a) with nodes labeled by Lorentzian Livine- 
Speziale coherent intertwiners. The Lorentzian coherent 
intertwiner is the product of Livine-Speziale coherent in- 
tcrtwiner with a Lorcntzian-geometry phase. 

For the approximation scheme, we consider the limit, 
where the Barbero-Immirzi parameter is taken to zero 
7 — > 0, and the spin of the boundary state is taken to in- 
finity j — > 00, keeping the size of the quantum geometry 
A ~ 7j finite and fixed. This limit corresponds to ne- 
glecting Planck scale discreteness and twisting effects, at 
large finite distances. It is interesting to notice that the 
same limit was considered in [51] in the context of loop 
quantum cosmology. In this limit, the two-point func- 
tion we obtain exactly matches the one obtained from 
Lorentzian Regge calculus [3"8] . 

Deriving the LQG correlation function at the level of 
a single spin foam vertex is certainly only a first step. 
Within the setting of a vertex expansion, an analysis of 



2 The finite and fixed area A corresponds to the finite and fixed 
distance between the two points where the correlation functions 
are defined. 
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the LQG correlation function for an arbitrary number of 
spinfoam vertices is needed. In [10] , asymptotic anal- 
ysis of spinfoams with an arbitrary number of vertices 
is studied in 7 — > limit. Without taking 7 — !> 0, the 
large- j of spinfoams with an arbitrary number of vertices 
is studied in [5] with a closed manifold, and is generalized 
recently to the one with boundary [TTJ [T2] . It would be 
interesting to investigate the contribution of the 7-term 
to correlation functions when more than a single spin 
foam vertex is considered. 
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